

function phi = vercheil(vec,l,dbete,deltaxx);
    
myfun = @(x) vercheil2(vec,l,dbete,x);
infun = @(x) integral(myfun,0,x)-deltaxx;

ttfun =@(x)myfun(x)*vec^2/l*tan(deg2rad(dbete*x))*sind(-vec/l/deg2rad(dbete)*log(cos(deg2rad(dbete*x))));
phi = newiteration(infun,ttfun,0.1,0.01);

end